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Nuclear forces can be systematically derived using effective chiral Lagrangians consis- 
tent with the symmetries of QCD. I review the status of the calculations for two- and 
three-nucleon forces and their applications in few-nucleon systems. I also address issues 



Qs ' like the quark mass dependence of the nuclear forces and resonance saturation for four- 



nucleon operators. 

1. INTRODUCTION 

The forces between two (a few) nucleons is one of the most studied problems in theo- 
retical physics. Over many decades, a standard picture had evolved, in which these forces 
are described in terms of meson exchanges (which is an extension of the range expansion 
suggested in the 1950ties by Taketani and collaborators [[I])- Nuclei, which comprise most 
^ ■ of the baryonic matter in the universe, are made of slowly moving nucleons, so that given 
a potential, one simply has to solve the Schrodinger equation for the A-nucleon system, 

H^ A = ElpA (1) 

to obtain the properties of nuclei. More precisely, the nuclear Hamiltonian can be written 

as 

H=T+V, V = J2 V v + V ^ + --- , (2) 

with T the kinetic energy operator and the potential V is a string of terms comprising 
two- and three-nucleon forces. Given the two-nucleon potential Vij obtained e.g. from the 
analysis of the many NN scattering data and adjusting a few parameters in the three- 
nucleon potential Vijk, the spectra for nuclei up to A ~ 12 (and other properties) can 
be calculated with high precision using Monte-Carlo methods leading to an astonishing 
agreement with data (see e.g. [0 El)- Other precision methods exist for few-nucleon 
systems, like e.g. exact solutions of the Faddeev-Yakubovsky equations, the stochastic 
variational method or the use of hyperspherical harmonics. However, in such a framework 
based on meson exchange or semi-phenomenological potentials, one can not explain why 
the 2N forces are so much stronger than the 3N ones and also, when including external 
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sources, gauge invariance is not easy to obey (but there are prescriptions to do so). 
Furthermore, the connection to the fundamental theory of the strong interactions, QCD, 
is loose and it is difficult to estimate the theoretical errors. On the other hand, over 
the last few decades Chiral Perturbation Theory (CHPT) has become a standard tool 
for analyzing the properties of hadronic systems at low energy where the perturbative 
expansion of QCD in powers of the coupling constant cannot be applied. It is based 
on the approximate and spontaneously broken chiral symmetry of QCD. Starting from 
the most general effective Lagrangian for Goldstone bosons (pions in the two-flavor case 
of u and d quarks) and matter fields (nucleons, . . .) consistent with the symmetries of 
QCD, the hadronic S-matrix elements are obtained via a simultaneous expansion in small 
external momenta and quark masses (where small refers to the scale of symmetry breaking, 
A x ~ 1 GeV) . Goldstone boson loops are incorporated to obey perturbative unitarity and 
all corresponding ultraviolet divergences can be absorbed at a given order in the chiral 
expansion by the counterterms of the effective Lagrangian. This perturbative scheme 
works well in the pion and the pion-nucleon sectors, where the interaction vanishes when 
the external momenta go to zero (in the chiral limit). In the case of a few interacting 
nucleons, the situation is very different in that one has to deal with a non-perturbative 
system. Indeed, perturbation theory is expected to fail already at low energy due to the 
presence of the shallow few-nucleon bound states. To make this more precise, observe 
that the np scattering length in the 1 Sq partial wave is much bigger than the largest 
natural scale set by the pion Compton wavelength, |a( 1 5'o)| — 24 fm ^> 1/M n ~ 1.4 fm 
or consider the binding momentum in the deuteron, which is much smaller than the pion 
mass, 7 = a/ Bdrrid/2 ~ 45 MeV <C M n ~ 140 MeV. A suitable non-perturbative approach 
has been suggested by Weinberg [Ej, who showed that the strong enhancement of the few- 
nucleon scattering amplitude arises from purely nucleonic intermediate states. Weinberg 
suggested to apply power counting to the kernel of the corresponding scattering equation, 
which can be viewed as an effective nuclear potential. This idea has been explored in 
the last decade by many authors. In the following I will show how far this program 
has matured in the description of few-nucleon systems. Space does not allow for giving 
an accurate representation of the historical developments and also I will not discuss the 
pionless effective field theory (EFT) here (for comprehensive reviews see [0E])- 

2. EFFECTIVE CHIRAL LAGRANGIAN 

Here I briefly discuss the effective chiral Lagrangian £ e g underlying the EFT calculation 
of the nuclear forces. The QCD Lagrangian for massless up and down quarks is invari- 
ant under global flavor SU(2)l x SU(2)r transformations or, equivalently, under vector 
and axial-vector transformations. This is called chiral symmetry. Among other facts, 
the absence of parity doublets of low mass hadrons suggests that the axial symmetry is 
spontaneously broken. The pions are the natural candidates for the required Goldstone 
bosons. They acquire a non-vanishing mass due to the explicit symmetry breaking caused 
by the small up and down quark masses. We are interested in low-energy nuclear physics 
here, where the degrees of freedom are the composite hadrons. Their interaction can de- 
scribed by an effective Lagrangian, which has to be constrained by chiral symmetry and 
should include explicitly symmetry breaking parts proportional to powers of the quark 
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masses (for the relation between the effective chiral Lagrangian and QCD, see [0IE])- For 
the applications to be considered, nucleon momenta are comparable to the pion mass and 
somewhat higher, but still smaller than the p-mass. In that case a standard one-boson 
exchange picture turns into NN contact forces for the heavy meson exchanges and only 
the few-pion exchanges are kept explicitly. The construction of the most general effective 
Lagrangian out of pion and nucleon fields constrained by chiral symmetry is by now stan- 
dard, based on the non-linear realization of chiral symmetry. There is an infinite number 
of possible terms, which can be ordered according to the parameter 

A = d+^n-2 (3) 

characterizing the vertices (note this differs from the standard counting done in the pion 
and pion-nucleon sectors). Here d is the number of derivatives and n the number of 
nucleon field operators. Spontaneously broken chiral symmetry enforces A > 0. The first 
few terms for the interacting effective Lagrangian after a p/m expansion take the form 
(for an early detailed review on the construction principles and applications, see [|5]) 



TCT ■ • V7T H — T ■ (TT X TT) + ■ ■ 



-eff ' 2F ■ - — ' AF 2 



N 



--C s (N*N) (N^N) - -C T (N^ctN) (rftrN) , 

^eff = +~E^ Nl [-2CiM^7T 2 + Czd^Txd^TT - \c A Ei jk E abc aiT a { V jlX b ) (Vfc7T c ) + 



N (4) 



° (N*N)(N*tr ■ tN) ■ -Vtt - -E (N*N)(N*tN) 2 , (5) 



4iV ' 2 

4? = ~C X [{N^VNf + (VMiV) 2 ] + ■ ■ ■ + C^N^aAN^^N) + ■ • • , (6) 

where the upper index refers to A = 0, 1 and 2. Here, N and tt denote the isodoublet 
nucleon and the isovector pion fields, in order. The parameters of £ e ff, the so called 
low-energy constants (LECs), can be partitioned in several groups: some can be de- 
termined in the n/n-N system (g^, F^, ci, C3, C4) and others from nucleonic systems only 
(Cs,Ct, E,Cx, ...C7). The constant D also affects the NN7r system, see [CD]. All these 
constants are of course not determined by chiral symmetry, but have to be adjusted to 
experimental data. 



3. THE HIERARCHY OF NUCLEAR FORCES 

The strength of the EFT approach is based on its underlying power counting, which 
allows to organize all possible contributions to the iV-nucleon potential in a systematic 
way. As a consequence of that, one naturally obtains a hierarchy of nuclear forces, as was 
shown already by Weinberg [0] and van Kolck [E]. Here, I will discuss how this comes 
about including corrections at next-to-next-to- next-to-leading order (N 3 LO). According 
to chiral power counting the dominant contributions to the effective Hamiltonian for few 
nucleons are of the order (Q/A)°, where Q ~ M„ refers to the soft scale (typical momenta 
involved in the process) and A to the hard scale (the chiral symmetry breaking scale, or 
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Figure 1. The hierarchy of the 
nuclear forces in chiral effective 
field theory. Shown are represen- 
tative diagrams (contact interac- 
tions, one-, two- and three-pion 
exchanges) symbolizing the var- 
ious topologies at leading, next- 
to-leading, . . . order. Due to par- 
ity, 2N contact interactions can 
only appear at even orders in 
Q/ A. The two-pion exchange di- 
agrams at N 2 LO (N 3 LO) include 
couplings with two (three) deriva- 
tives (or pion mass insertions) 
from the ttN effective chiral La- 
grangian. At N 3 LO, only a few 
topologies are shown as indicated 
by the ellipsis. 





2N forces 


3N forces 


4N forces 




X 









NLO {%) 


) 

/ 


( 


















N 2 LO {%) 














\ 


— ^ 

< 


i — 

X 




N 3 LO (£) 


) 


(1 


• 


• 


• 








+ 


• • • 


+ 

■ ■ • 



the ultraviolet cut-off to render the scattering equation finite). These leading-order (LO) 
contributions turn out to be of the 2N type, given by one-pion exchange (OPE) and two 
four-nucleon contact interactions without derivatives, cf. Fig. 1. The first corrections 
at next-to-leading order (NLO) are of the order (Q/A) 2 and still of the 2N type. They 
result from two-pion exchange (TPE) with the leading (i.e. with one derivative) nNN and 
ititNN vertices and seven independent NN contact interactions with two derivatives. The 
LECs entering the expressions for TPE at NLO are the nucleon axial-vector coupling 
and the pion decay constant F n . Both LECs are measured rather accurately, so that the 
leading TPE contribution is parameter-free. On the contrary, the LECs accompanying 
the four-nucleon contact operators are unknown and have to be fixed from a fit to low 
angular momentum partial waves. Thus, at LO and NLO we have only two-nucleon forces 
and therefore expect more-nucleon forces to be parametrically suppressed, in agreement 
with empirical information. 

At NNLO (~ (Q/A) 3 ) one has to take into account the subleading TPE contributions 
given by the triangle diagram with the 7T7tNN vertex with two derivatives or one insertion 
of M 2 . The corresponding LECs are denoted and have been fixed in the nN system, 
see e.g. [El EH E]- The numerical values of these LECs found in several analyses of irN 
scattering in CHPT are rather large compared to what is expected on dimensional reasons. 
Similar large values of 03,4 have also been obtained recently from the np and pp partial 
wave analysis carried out by the Nijmegen group [EI!- The large values of the 03,4 can at 
least be partially explained by the fact that these LECs are saturated by the A-excitation 
[IT2]. Other important contributions to C\ and C4 can be traced back to meson-exchange 
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in the t-channel. The large numerical values of the Cj's gives rise to a subleading TPE 
contribution to the NN potential that shows an unphysically strong attraction already at 
intermediate distances r ~ 1 — 2 fm when standard dimensional regularization in the pion 
loop integrals is applied. This was already pointed out in Ref. [QS]. To circumvent this 
problem, spectral function regularization was proposed in [QZj, as explained in the next 
section. It is important to note that at this order the first non-vanishing three-nucleon 
forces (3NF) appear, given by the three topologies shown in Fig. 1. I will come back to 
these later. Finally, at N 3 LO, one has to consider four-nucleon terms with four derivatives 
(or pion mass insertions), further corrections to the two-pion exchange (including now 
some of the dimension three couplings di from ) as well as the leading order three-pion 
exchange. These TPE and 3PE corrections were worked out explicitely utilizing DR in 
[ I18j. For the NN problem, one has 15 independent four-nucleon operators that feed 
into the S-, P- and D-waves and the corresponding mixing angles. Further corrections to 
the 3NF also appear at this order, it is however very important to stress that these are 
free of unknown six-nucleon LECs. In addition, the first non-vanishing corrections to the 
four-nucleon force appear at this order. Consequently, assuming that all coefficients are 
of natural size, chiral symmetry applied to the few-nucleon potential gives the hierarchy 
of the forces 

{V 2N ) A > (V m ) A > (Vw)a , (7) 
which is in agreement with phenomenological determinations of these forces. 

4. SPECTRAL FUNCTION REGULARIZATION 

Before discussing the physics related to few-nucleon systems, there is one important 
technical development to be discussed in a bit of detail here. As already stated, the 
numerically large values of the LECs q found in the tcN system lead to an unphysically 
strong attraction of subleading TPE. This is related to the use of dimensional regulariza- 
tion in the pion loop integrals. This problem can be overcome employing spectral function 
regularization, which has been used early in the construction of two-pion exchange con- 
tributions to the NN potential based on dispersion theory [d]. To be specific, consider 
the isoscalar central part of the subleading TPE which results from the triangle diagram 
and is given by 

Vc (q) = -M- / -^L , P ~ q2 _ (% Cl Ml + c 3 (/ 2 - q 2 )) , (8) 

where q is the nucleon momentum transfer and q = \q\, I = The integral is cubically 
divergent and needs to be regularized. Applying dimensional regularization (DR) one 
finds: 

Vc(q) = (2M 2 (2 Cl - c 3 ) - c 3 g 2 ) (2M 2 + g 2 )l arctan + . . . . ( 9 ) 

The ellipses refer to polynomial terms of the kind a + (3q 2 , whose explicit form is not of 
relevance here. In order to obtain the potential in coordinate space one has to make an 
inverse Fourier-transform of Vc(q) in Eq. Q. The ordinary inverse Fourier-transform is 
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Figure 2. The potential Vc in r-space. The solid line (shaded band) shows the DR 
(spectral function regularized, A = 500 . . .800 MeV) result. The dashed (dashed-dotted) 
line refers to the phenomenological a (cr + cj + p) contributions based on the isospin triplet 
configuration space version (OBEPR) of the Bonn potential. 



obviously not possible due to the the fact that Vc(q) grows with q. One can nevertheless 
obtain Vc(r) at each r > using the spectral function representation [Ho]: 



V c (q) = — f dp\ 



p 2 + q 2 



(10) 



where the spectral function can be obtained from Vc(q) in Eq. Q via 



p{p) = 3 [V c (0 + - ip 



^ (2M^(2 Cl - c 3 ) + c 3 p 2 ) (2M 2 - /i 2 )^(/i - 2M W ) . (11) 



In eq. ffTUj) a twice subtracted dispersion integral is given which is needed in order to 
account for the large-/i behavior of p(/i). 

The inverse Fourier-transform in terms of the spectral function p(p) can easily be evalu- 
ated via 



Vc(r) 



2n 2 r 



dppe Mr p(/x). 



(12) 



2M n 



Substituting p(p) from Eq. (fTT|) into Eq. (fT2"j) and using for the LECs ci^ the central 
values from [Ej, c\ = —0.81 GeV -1 and c 3 = —4.70 GeV -1 , one obtains the coordinate 
space representation of the potential Vc(r) shown by the solid line in Fig. |21 The central 
part of the NNLO TPE potential turns out to be several times more attractive at inter- 
mediate distances than the phenomenological a (a + u + p) contributions. This unphysical 
attraction shows up in the D and F-wave phase shifts as observed in [I16|. The origin 
of the unphysical attraction at NNLO can be traced back by looking at the integral in 
Eq. f[T2"j) . While at large distances the integral is dominated by low-/i components (of 
the order p ~ 350 MeV), already at intermediate distances rather high-/i terms (of the 
order p ~ 600 MeV) in the spectral function provide a dominant contribution. Clearly, at 
shorter distances even higher-/i components become important. Chiral EFT can hardly 
provide convergent results for the spectral function at p ~ 600 MeV and higher. Instead 
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of keeping such large-// contributions in the regularized loop integral expressions it is ad- 
vantageous to perform the integration in Eqs. (|1UI12|) only over the low-// region, where 
chiral EFT is applicable. This can be achieved by introducing the regularized spectral 
function 

p(/i)^p A (//) = p(/i)fl(A-/x), (13) 

with the reasonably chosen cut-off A < M p . Certainly, taking a too small A in Eq. ()13|) will 
remove the truly long-distance physics while too large values for the cut-off may affect the 
convergence of the EFT expansion due to the inclusion of spurious short-distance physics. 
In Fig. |21 we show Vc(r) obtained using the spectral function regularization Eq. ()13|) with 
A = 500 . . . 800 MeV. The strongest effects of the cut-off are observed at intermediate 
and short distances, where the unphysical attraction in dimensionally regularized TPE 
is greatly reduced. On the other hand, the asymptotic behavior of the potential at large 
r is not affected by the choice of regularization. The results for D- and F-waves are 
greatly improved when the spectral function regularization is used instead of DR, see [ 
IT7] . One can also show that the spectral function regularization is equivalent to (finite) 
momentum cut-off regularization of pion loop integrals. It should be understood that 
this new regularization scheme does not introduce any model dependence in the EFT 
procedure as soon as A is chosen of the order of (or larger than) M p . Various choices for 
A (including A = oo, which is equivalent to DR) differ from each other by higher-order 
contact terms and lead to exactly the same result for observables provided one keeps terms 
in all orders in the EFT expansion. Of course, this choice of regularization in Eq. (fT3~j) is 
by no means unique. Different choices lead to equivalent results for the potential up to 
higher order terms and may be used as well. The advantage of the form Eq. (J13|) is that 
it does not generate spurious long-range contributions which are suppressed by inverse 
powers of A. For a more detailed discussion on this issue, I refer to [ITTj. 

5. TWO NUCLEONS AT N 3 LO 

The interactions between two nucleons at N 3 LO were studied in [I2D1- Here, I can 
only sketch the most salient features of that study. As discussed before, at this order 
the two-nucleon potential consists of one-, two- and three-pion exchanges and a set of 
contact interactions with zero, two and four derivatives, respectively. We have applied 
spectral function regularization to the multi-pion exchange contributions. Within this 
framework, we have shown that three-pion exchange can safely be neglected. The corre- 
sponding cut-off is varied from 500 to 700 MeV. The LECs related to the dimension two 
and three NNtttt vertices are taken consistently from studies of pion-nucleon scattering 
in chiral perturbation theory [E3 [T^j. in the isospin limit, there are 24 LECs related 
to four-nucleon interactions which feed into the S-, P- and D-waves and various mixing 
parameters. In addition, various isospin breaking mechanisms were considered. In the ac- 
tual calculations, we have included the leading charge- independence and charge-symmetry 
breaking four-nucleon operators, the pion and nucleon mass differences in the 1PE, and 
the same electromagnetic corrections as done by the Nijmegen group (the static Coulomb 
potential and various corrections to it, magnetic moment interactions and vacuum po- 
larization). This is done because we fit to the Nijmegen partial waves. In the future, 
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it would be important to also include mass differences in the 2PE, ^-exchange and the 
isospin breaking corrections to the pion-nucleon scattering amplitude (which have been 
consistently determined in [I21j). We therefore have phases for the pp, np and nn systems. 
To obtain the bound and scattering states, we use the Lippmann-Schwinger equation with 
the relativistic form of the kinetic energy. Such an approach can easily be extended to 
external probes or few-nucleon systems. The LS equation is regulated in the standard 
way, namely by multiplying the potential V(p, p') with a regulator function f A , 

V(p, p') -> f A (p) V(p, f) f\p') . (14) 

We use the exponential regulator function f A (p') = exp[— p 6 /A 6 ] , with the cut-off varied 
from 450 to 600 MeV. The total of 26 four-nucleon LECs has been determined by a 
combined fit to some np and pp phase shifts from the Nijmegen analysis together with 
the nn scattering length value a nn = — 18.9fm. The resulting LECs are of natural size 
except D\ SQ and D^ S1 . In contrast to the fits at NLO and NNLO, we had to extend 
the fit range to E\ ah = 200 MeV. The description of the low phase shifts (S, P, D) is 
excellent, see Fig. El for the S- and some selected P- and D- waves. In all cases, the N 3 LO 




Lab. Energy [MeV] Lab. Energy [MeV] Lab. Energy [MeV] 

Figure 3. Selected np phase shifts and mixing angle €2 versus the nucleon laboratory 
energy. The grid, light shaded and dark shaded bands show the NLO, NNLO [ IT7| and 
N 3 LO results, respectively. The filled circles depict the Nijmegen PWA results [|22j and 
the open triangles are the results from the Virginia Tech PWA [ I23j . 



result is better than the NNLO one with a sizeably reduced theoretical uncertainty. This 
holds in particular for the problematic 3 Pq wave which was not well reproduced at NNLO. 
The peripheral waves (F, G, H, . . .), that are free of parameters, are also well described 
with the expected theoretical uncertainty related to the cut-off variations, see Fig. |3] for 
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e [deg] 6 [deg] 

Figure 4. np differential cross section (upper row) and vector analyzing power (lower 
row) at Eiab = 50 MeV (left panel) and = 96 MeV (right panel). Shaded bands refer 
to the NNLO result, dashed lines to the Nijmegen phase shift analysis (NPSA) [ 22 . 



1 F 3 and l G±. We stress that the description of the phases in general improves when 
going from LO to NLO to NNLO to N 3 LO, as it is expected in a converging EFT. The 
resulting S-wave scattering lengths and range parameters in the np and pp systems are 
in good agreement with the ones obtained in the Nijmegen PWA. In addition, we can 
give theoretical uncertainties for all these quantities, which are mostly in the one percent 
range. The scattering observables (differential cross sections, analyzing powers) for the np 
system as shown e.g. in Fig. 0] are well described, with a small theoretical uncertainty at 
the order considered here. The deuteron properties are further predictions. In particular, 
we have not included the binding energy in the fits, the deviation from the experimental 
value is in the range from 0.4 to 0.07%. The asymptotic S-wave normalization and the 
asymptotic D/S ratio are also well described, see Table 1. The remaining discrepancies 
in the quadrupole moment and the rms matter radius are related to the short-ranged 
two-nucleon current not considered in [I2D]- Note that the 2N system at this order was 
also studied in [121], utilizing dimensional regularization in the pion loop graphs. In that 
work, however, no detailed dicsussion of the theoretical uncertainties was given. 



Table 1 

Deuteron properties at various orders in chiral EFT in comparison to the data. 





NLO 


N 2 LO 


N 3 LO 


Exp. 


E d [MeV] 


-2.171 ...-2. 188 


-2.189 ... -2.202 


-2.216 ...-2.223 


-2.2246 


A s [im 1 ' 2 } 


0.868... 0.873 


0.874... 0.879 


0.882. . .0.883 


0.8846(9) 


V 


0.0256 ...0.0257 


0.0255 ...0.0256 


0.0254 ...0.0255 


0.0256(4) 
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6. THREE NUCLEONS AT N 2 LO 

In the chiral EFT framework, three-nucleon forces (3NFs) arise consistently with the 
iViV forces from the effective Lagrangian. This is one of the major advantages of this ap- 
proach - such a consistency has never been achieved before. Furthermore, chiral EFT also 
allows to generate the most general structures consistent with the underlying symmetries, 
this is simply a consequence of utilizing the most general effective Lagrangian consistent 
with these principles. As noted earlier, the leading non-vanishing contributions to the 
3NF appear at NNLO, given by three different topologies (see Fig. 1): the TPE graphs, 
the OPE diagram and the six-nucleon contact interactions. The TPE contribution is 
given in terms of LECs 01^4 from the pion-nucleon system (see [ I25| I26|) 

T/3NF _ V- 1 ( 9A\ 2 (^i-Qi)(^j-Qj) F a0a0 ,x 

^ "2^2 \2Fj (q^ + M^ + M^^ T * ' 

where Qi = Pi ' — Pi, Pi {Pi ') are the initial (final) momenta of the nucleon % and 



4 Cl 7l4 2 2c 3 

7T J 7T 
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The OPE and contact contributions are given in terms of the LECs D and E, respectively, 
and take the form 



y3NF 
cont 



WE{r r r k ). (17) 



2 

3+k 



To pin down the 3NF at this order, one needs two observables to determine the LECs E 
and D. In these two parameters were determined from the 3 H binding energy and 
the 2 a n d scattering length. In that paper, a detailed study of 3- and 4-nucleon observables 
including this 3NF was performed. The role of the 3NF increases with increasing energy, as 
exemplified in Fig. 5, where the minimum of the differential cross section for nd scattering 
at 65 MeV is shown. It is clearly visible that the calculation without 3NF badly misses 
the data. I should stress again that the solid curve including the 3NF is a prediction since 
all parameters have been fixed before. 

In Fig. IH1 the tensor analyzing powers T 20 and T 2 2 for elastic pd scattering at 70 MeV are 
shown [I23I2H1- While T 20 is well described by the NNLO calculation, one observes some 
visible differences in T 22 . It is conceivable that this will be cured by the N 3 LO 3NFs, which 
are presently worked out. We believe that these corrections to the 3NF are free of LECs, 
that means even at N 3 LO one has only two parameters to completely pin down the chiral 
3NF. There is also a large amount of pd break-up data, for a detailed comparison between 
theory and experiment, see [EE]- For illustration, in Fig. the break-up cross section 
in the symmetric space-star configuration as well as the analyzing power are shown, in 
comparison to the conventional approach based on two high-precision iViV potentials and 
two available phenomenological 3NFs. In [l2lfl. we also solved the Yakubowsky equations 
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Figure 5. Differential cross sec- 
tion for nd scattering at 65 MeV 
around the minimum at 9 ~ 130°. 
Solid (dashed) line: NNLO pre- 
diction with (without) 3NF. The 
dot-dashed line gives the predic- 
tion based on the CD-Bonn po- 
tential without 3NF. 




60 120 180 60 120 180 



Bern, [deg] c . m . [deg] 

Figure 6. Left: Tensor analyzing powers T 20 (left) and T 2 2 (right) in pd scattering at 
70 MeV recently measured at RIKEN [I2E|- The green band shows the NNLO calculation. 

to determine the binding energy (BE) of the a-particle. One finds (all numbers in MeV): 

NNLO: BE( 4 He) = -29.98. ..- 29.51 , Exp.: BE( 4 He) = -29.8 ± 0.1 , (18) 

where the experimental value is the "synthetic" binding energy for pure np forces to allow 
for a direct comparison with the then available chiral EFT forces (for details see [ I26j). 
There has also been some more recent work on the isospin dependence of the three-nucleon 
force and its effect on the binding in the 3-nucleon system [ EIH EI] ■ 

7. RESONANCE SATURATION OF FOUR-NUCLEON COUPLINGS 

In the meson sector of CHPT it has been demonstrated that the numerical values 
of the pertinent NLO LECs Li can be understood to a high precision by integrating 
out heavy mesons of all types from the theory [ ESj • This was coined resonance 
saturation. For the dimension two and and some dimension three LECs of the pion- 
nucleon sector, these ideas were extended in [ |33] and later used in studies of neutral 
pion electroproducion off protons. It thus appears natural to confront the four-nucleon 
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s [MeV] S [MeV] 

Figure 7. Left: pd breakup cross section data in [mb MeV -1 sr -2 ] along the kinematical 
locus S (in MeV) at 65 MeV. The data are from NNLO predictions (dark shaded 

band) and NLO (light shaded band) compared to the conventional NN forces+TM 3NF 
predictions. The solid (dashed) line refers to the AV18+URBANA IX (CD Bonn+TM') 
results. Right: Analyzing power. 



LECs determined from nuclear EFT with the highly successful phenomenological/meson 
exchange models of the nuclear force following the lines of Ref. [ EH] and extending the 
ideas of resonance saturation to this sector. To be specific, consider some genuine one- 
boson-exchange (OBE) models. In these models the long range part of the interaction is 
given by OPE (including a pion-nucleon form factor) whereas shorter distance physics is 
expressed as a sum over heavier meson-exchange contributions: 

vnn = k+ Yl Vm - ( 19 ) 

M=a,p,... 

Here some mesons can be linked to real resonances (like e.g. the p-meson) or are param- 
eterizations of certain physical effects, e.g. the light scalar-isoscalar a-meson is needed 
to supply the intermediate range attraction (but it is not a resonance). The correspond- 
ing meson-nucleon vertices are given in terms of one (or two) coupling constant (s) and 
corresponding form factor (s), characterized by some cut-off scale. These form factors are 
needed to regularize the potential at small distances (large momenta) but they should not 
be given a physical interpretation. As depicted in Figd (left panel) for nucleon momen- 
tum transfer below the masses of the exchanged mesons, one can interpret such exchange 
diagrams as a sum of local operators with increasing number of derivatives (momentum 
insertions). This is explained in detail in Ref. [EH]- in that work the short-range part of 
different phenomenological potential models was power expanded and the resulting con- 
tact operators were compared with the ones in the EFT approach. The latter have to be 
corrected by adding the corresponding power expanded TPE contributions, which are not 
present in most phenomenological models. It was then demonstrated explicitly that the 
values of the LECs Cj determined from various phenomenological OBE models are close 
to the values found in EFT at NLO and NNLO, see Fig. [7| This was repeated in [ IT7] 
utilizing spectral function regularization with a similar outcome - the LECs appearing in 
the chiral EFT fits can be well represented by heavy meson exchanges. 
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Figure 8: Top: Heavy meson ex- 
change represented as a sun of 
local interactions with increasing 
number of derivatives. Right: 
Comparison of the NLO (green 
bands) and NNLO (red bands) 
LECs with values obtained from 
accurate potentials [I3E]. 
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8. QUARK MASS DEPENDENCE OF THE NUCLEAR FORCES 

Because of the smallness of the up and down quark masses, one does not expect sig- 
nificant changes in systems of pions or pions and one nucleon when the quark masses 
are set to zero (with the exception of well understood chiral singularities like e.g. in 
the pion radius or the nucleon polarizabilities). The situation is more complicated for 
systems of two (or more) nucleons, as first discussed in EFT in [ EE] . Here, I report on 
similar work [ EI] that is mostly concerned with the properties of the deuteron and the 
S-wave scattering lengths as a function of the quark (pion) mass. These questions are not 
only of academic interest, but also of practical use for interpolating results from lattice 
gauge theory. E.g. the S-wave scattering lengths have been calculated on the lattice using 
the quenched approximation [EH]. Another interesting application is related to imposing 
bounds on the time-dependence of some fundamental coupling constants from the NN 
sector, as discussed in [ESJ. To address this issue, at NLO the following contributions 
have to be accounted for (in addition to the LO OPE and contact terms without deriva- 
tives): i) contact terms with two derivatives or one M^-insertion, ii) renormalization of 
the OPE, iii) renormalization of the contact terms, and iv) two-pion exchange (TPE). 
This induces explicit and implicit quark mass dependences. In the first category are the 
pion propagator that becomes Coulomb-like in the chiral limit or the M% corrections to 
the leading contact terms. These are parameterized by the LECs Ds,t at NLO. These 
LECs can at present only be estimated using dimensional analysis and resonance satura- 
tion [ 35J. The implicit pion mass dependence enters at NLO through the pion-nucleon 
coupling constant (note that the quark mass dependence of the nucleon mass only enters 
at NNLO) expressed through the pion mass dependence of gA/F n in terms of the quantity 

a ={Sfi - t> + (M - - ^ - SS 111 1 ' (20) 

Here U, d\s and die are LECs related to pion and pion-nucleon interactions, and the value 
of the pion mass is denoted by in order to distinguish it from the physical one denoted 
by M n . In particular, d 16 has been determined in various fits to describe nN — ► imN 
data, see [EOj. The deuteron BE as a function of the pion mass is shown in Fig El we find 
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that the deuteron is stronger bound in the chiral limit (CL) than in the real world, 

Bg L = 9.6 ± 1.9^ ' MeV , (21) 

where the the first indicated error refers to the uncertainty in the value of Dss 1 and 
di e being set to its average value while the second indicated error shows the additional 
uncertainty due to the uncertainty in the determination of di§ We find no other bound 
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Figure 9. Left panel: Deuteron BE versus the pion mass. The shaded areas show allowed 
values. The light shaded band corresponds to our main result with the uncertainty due 
to the unknown LECs Ds s t- The dark shaded band gives the additional uncertainty due 
to the uncertainty of d^. The heavy dot shows the BE for the physical case M n = M v 
Right panel: The inverse S-wave scattering lengths as functions of M n . The shaded areas 
represent the allowed values according to our analysis. The heavy dots corresponds to 
the values in the real world. The triangles refer to lattice QCD results from 



states, although the higher S = 1 partial waves rise linear with momentum due to the 
Coulomb-like pion propagator. Last but not least, we found smaller (in magnitude) and 
more natural values for the two S-wave scattering lengths in the chiral limit, 

a C L( 1 5 ) = -4.1±1.6^-°fm, and a CL ( 3 Si) = 1.5 ± 0.4^ fm . (22) 

As stressed in [ |37] , one needs lattice data for pion masses below 300 MeV to perform a 
stable interpolation to the physical value of M n , cf the right panel in Fig. O We conclude 
that nuclear physics in the chiral limit is much more natural than in the real world. 
Another interesting application of the quark mass of the nuclear forces is the recently 
conjectured infrared renormalization group limit cycle in the three-nucleon system [I41j. 

9. BRIEF SUMMARY AND OUTLOOK 

From the previous sections it should have become obvious that few-nucleon systems can 
be studied in chiral effective field theory in a systematic and model-independent way. The 
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Figure 10: 2N and 3iV cur- 
rents. First row: Leading 
pion exchange graphs. Sec- 
ond row: Corrections of var- 
ious ranges. Third row: 3iV 
currents. Solid, dashed, and 
wiggly lines represent nucle- 
ons, pions and photons, in or- 
der. Only some representa- 
tive diagrams are depicted. 
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two-nucleon system has been analyzed at N 3 LO and accurate results for the deuteron and 
scattering observables have been obtained. Furthermore, 3N, 4N and even 6N systems 
[ 1321 have been studied at N 2 LO. For the first time, the chiral three-nucleon force has 
been consistently included and the results to this order look very promising. Many other 
applications have also been performed (some times using the so-called hybrid approach 
which utilizes the kernel from EFT and wave functions from semi-phenomenological po- 
tentials). I mention just a few here: pion-deuteron scattering [133]. neutral pion photo- 
and electroproduction on the deuteron 



133] [ E2| , Compton scattering off deuterium 
nuclear parity violation [ HJ] or solar fusion and the hep process [ HH] . Other interest- 
ing developments are related to the lattice, just to name a few examples of recent work: 
nuclear lattice simulations [EH], the two-nucleon potential in partially quenched lattice 
QCD [|5U] or the discussion of the size of the lattice required to simulate the two-nucleon 
system [1ST]. 

From my point of view, the following problems should be worked out next: 

• To study few-nucleon physics at N 3 LO, we have to work out the 3NF and 4NF at 
this order. It is gratifying to notice that no new six-nucleon contact interactions 
appear at this order and thus one has large predictive power. Furthermore, the 4NF 
first appears at this order and is parameter-free. It is expected that this force will 
be much smaller than the 3NF. Still, it will be interesting to obtain the quantitative 
size of this force. 

• Similarly, the currents corresponding to electroweak probes have to be constructed 
to the same order. The strength of the chiral EFT approach is the consistency 
with the forces and the automatic incorporation of gauge symmetry. Lots of work 
on this problem has already been done by the Korean group [ E2] (see also [ EH] 
54 ). However, most of these results need to be rederived using the method of 
unitary transformation to make them consistent with the forces. Some pertinent 
diagrams are shown in Fig. 10. Furthermore, three-nucleon currents have never 
been considered. 

The methods described here pave the way to a precision nuclear physics consistent with 
the symmetries of QCD and thus offer a sound theoretical foundation for this decade old 
problem. 
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